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Abstract 

Capacity scaling laws are analyzed in an underwater acoustic network with n regularly located nodes on a 
square, in which both bandwidth and received signal power can be limited significantly. A narrow-band model is 
assumed where the carrier frequency is allowed to scale as a function of n. In the network, we characterize an 
attenuation parameter that depends on the frequency scaling as well as the transmission distance. Cut-set upper 
bounds on the throughput scaling are then derived in both extended and dense networks having unit node density and 
unit area, respectively. It is first analyzed that under extended networks, the upper bound is inversely proportional 
to the attenuation parameter, thus resulting in a highly power-limited network. Interestingly, it is seen that the 
upper bound for extended networks is intrinsically related to the attenuation parameter but not the spreading factor 
On the other hand, in dense networks, we show that there exists either a bandwidth or power limitation, or both, 
according to the path-loss attenuation regimes, thus yielding the upper bound that has three fundamentally different 
operating regimes. Furthermore, we describe an achievable scheme based on the simple nearest-neighbor multi-hop 
(MH) transmission. We show that under extended networks, the MH scheme is order-optimal for all the operating 
regimes. An achievability result is also presented in dense networks, where the operating regimes that guarantee 
the order optimality are identified. It thus turns out that frequency scaling is instrumental towards achieving the 
order optimality in the regimes. Finally, these scaling results are extended to a random network realization. As a 
result, vital information for fundamental limits of a variety of underwater network scenarios is provided by showing 
capacity scaling laws. 
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I. Introduction 

Gupta and Kumar's pioneering work [1] characterized the connection between the number n of nodes 
and the sum throughput in a large-scale wireless radio network. They showed that the total throughput 
scales as Q{^^/n/\ogn) when a multi-hop (MH) routing strategy is used for n source-destination (S-D) 
pairs randomly distributed in a unit areaj^MH schemes are then further developed and analyzed in [3]-[9], 
while their throughput per S-D pair scales far slower than 0(1). Recent results [10], [11] have shown 
that an almost linear throughput in the radio network, i.e. 0(n^~^) for an arbitrarily small e > 0, which 
is the best we can hope for, is achievable by using a hierarchical cooperation (HC) strategyH Besides the 
schemes in [10], [11], there have been other studies to improve the throughput of wireless radio networks 
up to a linear scaling in a variety of network scenarios by using novel techniques such as networks with 
node mobility [12], interference alignment [13], and infrastructure support [14]. 

Together with the studies in terrestrial radio networks, the interest in study of underwater networks has 
been growing, due to recent advances in acoustic communication technology [15]-[18]. In underwater 
acoustic communication systems, both bandwidth and received signal power are severely limited owing 
to the exponential (rather than polynomial) path-loss attenuation with long propagation distance and 
the frequency-dependent attenuation. This is a main feature that distinguishes underwater systems from 
wireless radio links. Hence, the system throughput is affected by not only the transmission distance but 
also the useful bandwidth. Based on these characteristics, network coding schemes [17], [19], [20] have 
been presented for underwater acoustic channels, while network coding showed better performance than 
MH routing in terms of reducing transmission power. MH networking has further been investigated in 
other simple but realistic network conditions that take into account the practical issues of coding and 
delay [21], [22]. 

One natural question is what are the fundamental capabilities of underwater networks in supporting a 
multiplicity of nodes that wish to communicate concurrently with each other, i.e., multiple S-D pairs, 
over an acoustic channel. To answer this question, the throughput scaling for underwater networks was 
first studied [23], where n nodes were arbitrarily located in a planar disk of unit area, as in [1], and 
the carrier frequency was set to a constant independent of n. That work showed an upper bound on the 
throughput of each node, based on the physical model [1], which scales as n-^/o'^-^oioin v^))^ where a 
corresponds to the spreading factor of the underwater channel and Wq represents the branch zero of the 
Lambert W function [24] H Since the spreading factor typically has values in the range 1 < a < 2 [23], 
the throughput per node decreases almost as 0{n^^^") for large enough n, which is considerably faster 
than the Q(y/n) scaling characterized for wireless radio settings [1]. 

In this paper, capacity scaling laws for underwater networks are analyzed in two fundamental but 
different networks: extended networks [4], [5], [10], [25], [26] of unit node density and dense networks [1], 
[6], [10] of unit area0 Unlike the work in [23], the information-theoretic notion of network capacity is 
adopted in terms of characterizing the model for successful transmission. Especially, we are interested in 
the case where the carrier frequency scales as a certain function of n in a narrow-band model. Such an 
assumption leads to a significant change in the scaling behavior owing to the attenuation characteristics. 
Recently, the optimal capacity scaling of wireless radio networks has been studied in [28], [29] according 
to operating regimes that are determined by the relationship between the carrier frequency and the number 
n of nodes. The frequency scaling scenario of our study essentially follows the same arguments as those 

'We use the following notation: i) fix) = 0{g{x)) means that there exist constants C and c such that f{x) < Cg{x) for all x > c. 
ii) fix) = o(g(x)) means that lim 44 = 0. iii) f(x) = nig{x)) if g{x) = 0(/(s)). iv) f(x) = u{g{x)) if g{x) = o{f{x)). v) 
fix) = eig{x)) if fix) = 0{g{x)) and g{x) = 0{f{x)) [2]. 

^Note that the HC scheme deals with a subtle issue around quantization, which is not our main concern in this work. 

^The Lambert W function is defined to be the inverse of the function z = W{z)e^''^'' and the branch satisfying > —1 is denoted 

by Wo{z). 

*Since the two networks represent both extreme network realizations, a realistic one would be in-between. In wireless radio networks, 
the work in [27] generalized the results of [10] to the case where the network area can scale polynomially with the number n of nodes. In 
underwater networks, we leave this issue for further study. 
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in [28], [29]. We aim to study both information-theoretic upper bounds and achievable rate scalings while 
allowing the frequency scaling with n. To the best of our knowledge, such an attempt has never been 
done before in underwater networks. 

We explicitly characterize an attenuation parameter that depends on the transmission distance and also 
on the carrier frequency, and then identify fundamental operating regimes depending on the parameter. For 
networks with n regularly distributed nodes, we derive upper bounds on the total throughput scaling using 
the cut-set bound. In extended networks, our upper bound is based on the characteristics of power-limited 
regimes shown in [10]. We show that the upper bound is inversely proportional to the attenuation parameter. 
This leads to a highly power-limited network for all the operating regimes, where power consumption is 
important in determining performance. Interestingly, it is seen that contrary to the case of wireless radio 
networks, our upper bound heavily depends on the attenuation parameter but not on the spreading factor 
(corresponding to the path-loss exponent in wireless networks). On the other hand, in dense networks, our 
upper bound basically follows the arguments, similarly as in [27]: there exists either a bandwidth or power 
limitation, or both, according to the operating regimes (i.e., path-loss attenuation regimes). Specifically, the 
network is bandwidth-limited as the path-loss attenuation is small. This is because network performance 
on the total throughput is roughly linear in the bandwidth. However, at the medium attenuation regime, the 
network is both bandwidth- and power-limited since the amount of available bandwidth and received signal 
power affects the performance. Finally, the network becomes power-limited as the attenuation parameter 
increases exponentially with respect to more than ^/n, i.e., at the high attenuation regime. Hence, our 
results indicate that the upper bound for dense networks has three fundamentally different operating 
regimes according to the attenuation parameter. In addition, to show constructively our achievability 
result for extended networks, we describe the conventional nearest-neighbor MH transmission [1], and 
analyze its achievable throughput. We show that under extended networks, the achievable rate scaling 
based on the MH routing exactly matches the upper bound on the capacity scaling for all the operating 
regimes. An achievability is also presented in dense networks by utilizing the existing MH routing scheme 
with a slight modification — we identify the operating regimes such that the optimal capacity scaling is 
guaranteed. Therefore, this key result indicates that frequency scaling is instrumental towards achieving 
the order optimality in the regimes. Furthermore, a random network scenario is discussed in this work. 
We show that under extended random networks, the conventional MH-based achievable scheme is not 
order-optimal for any operating regimes. 

The rest of this paper is organized as follows. Section |ll] describes our system and channel models. 
In Section [nil the cut-set upper bounds on the throughput are derived. In Section |IVl the achievable 
throughput scalings are analyzed. These results are extended to the random network case in Section W\ 
Finally, Section rvTl summarizes the paper with some concluding remarks. 

Throughout this paper the superscript H, [■]ki, and || • ||2 denote the conjugate transpose, the (fc,i)-th 
element, and the largest singular value, respectively, of a matrix. I„ is the identity matrix of size n x n, 
tr(-) is the trace, det(-) is the determinant, and \X\ is the cardinality of the set X. C is the field of 
complex numbers and E[-] is the expectation. Unless otherwise stated, all logarithms are assumed to be 
to the base 2. 

II. System and Channel Models 

We consider a two-dimensional underwater network that consists of n nodes on a square such that 
two neighboring nodes are 1 and units of distance apart from each other in extended and dense 

networks, respectively, i.e., a regular network [25], [26]. We randomly pick a matching of S-D pairs, so 
that each node is the destination of exactly one source. Each node has an average transmit power constraint 
P (constant), and we assume that the channel state information (CSI) is available at all receivers, but not 
at the transmitters. It is assumed that each node transmits at a rate T(n)/n, where T(n) denotes the total 
throughput of the network. 

Now let us tum to channel modeling. We assume frequency-flat channel of bandwidth W Hz around 
carrier frequency /, which satisfies / ^ W, i.e., narrow-band model. This is a highly simplified model. 
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but nonetheless one that suffices to demonstrate the fundamental mechanisms that govern capacity scal- 
ing. Assuming that all the nodes have perfectly directional transmissions, we also disregard multipath 
propagation, and simply focus on a line-of-sight channel between each pair of nodes used in [10], [11], 
[27]. An underwater acoustic channel is characterized by an attenuation that depends on both the distance 
Tki between nodes i and k {i,k e {I, - ■ ■ ,n}) and the signal frequency /, and is given by 

A{n,J) = cor^AfY'' (1) 

for some constant cq > independent of n, where a is the spreading factor and a(/) > 1 is the 
absorption coefficient [16]. For analytical tractability, we assume that the spreading factor a does not 
change throughout the network, i.e., that it is the same from short to long range transmissions, as in 
wireless radio networks [1], [4], [10]. The spreading factor describes the geometry of propagation and is 
typically 1 < a < 2 — its commonly used values are « = 1 for cylindrical spreading, a = 2 for spherical 
spreading, and a = 1.5 for the so-called practical spreading. Note that existing models of wireless networks 
typically correspond to the case for which a{f) = 1 (or a positive constant independent of n) and a > 20 
A common empirical model gives a(/) in dB/km for / in kHz as [16], [30]: 

101oga(/) = ao + ai/^ + 027— — ^ + asT— — (2) 

where {ao, ■ ■ • ,03, 61, 62} are some positive constants independent of n. As stated earlier, we will allow 
the carrier frequency / to scale with the number n of nodes. As a consequence, a wider range of both 
/ and n is covered, similarly as in [27]-[29]. In particular, we consider the case where the frequency 
scales at arbitrarily increasingrates relative to n, which enables us to really capture the dependence on 
the frequency in performance|j The absorption a(/) is then an increasing function of / such that 



a{f) = e (e^^^') (3) 



with respect to / for some constant ci > independent of n. 

The noise rii at node i G {1, ■ ■ ■ , n} in an acoustic channel can be modeled through four basic sources: 
turbulence, shipping, waves, and thermal noise [16]. We assume that rij is the circularly symmetric complex 
additive colored Gaussian noise with zero mean and power spectral density (PSD) N{f), and thus the 
noise is frequency-dependent. The overall PSD of four sources decays linearly on the logarithmic scale in 
the frequency region 100 Hz - 100 kHz, which is the operating region used by the majority of acoustic 
systems, and thus is approximately given by [16], [31] 

logiV(/) = a4-a5log/ (4) 

for some positive constants 04 and 05 independent of This means that N{f) = 0(1) since 

Nif) = e f (5) 



in terms of / increasing with n. From (|3]) and dS), we may then have the following relationship between 
the absorption a{f) and the noise PSD N{f): 

Nif) = (TTZZzh^) ■ (6) 



.(loga(/))«5/2. 

From the narrow-band assumption, the received signal yk at node k G {1, ■ ■ ■ , ra} at a given time 
instance is given by 

Vk = y^^hkiXj + nk, (7) 

^The counterpart of a in wireless radio channels is the path-loss exponent. 

^Otherwise, the attenuation parameter a(/) scales as 0(1) from (|2j, which is not a matter of interest in this work. 
^Note that in our operating frequencies, 05 — 1.8 is commonly used for the above approximation [16]. 
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where 



hki 



(8) 



VMrki, f) 



represents the complex channel between nodes i and k, Xi G C is the signal transmitted by node i, and 
/ C {I,-- - ,n} is the set of simultaneously transmitting nodes. The random phases 9ki are uniformly 
distributed over [0, 27i) and independent for different i, k, and time. We thus assume a narrow-band 
time-varying channel, whose gain changes to a new independent value for every symbol. Note that this 
random phase model is based on a far-field assumption [10], [11], [27]|!| which is valid if the wavelength 
is sufficiently smaller than the minimum node separation. 

Based on the above channel characteristics, operating regimes of the network are identified according 
to the following physical parameters: the absorption a(/) and the noise PSD A^(/) which are exploited 
here by choosing the frequency / based on the number n of nodes. In other words, if the relationship 
between / and n is specified, then a(/) and N(f) can be given by a certain scaling function of n from 
© and ©, respectively. 



To access the fundamental limits of an underwater network, new cut- set upper bounds on the total 
throughput scaling are analyzed from an information-theoretic perspective [33]. Consider a given cut L 
dividing the network area into two equal halves, as in [10], [27] (see Figs. \T\ and |2] for extended and 
dense networks, respectively). Under the cut L, source nodes are on the left, while all nodes on the right 
are destinations. In this case, we have an Q{n) x Q(n) multiple-input multiple-output (MEMO) channel 
between the two sets of nodes separated by the cut. 

A. Extended Networks 

In this subsection, an upper bound based on the power transfer argument [10] is established for extended 
networks, where the information flow for a given cut L is proportional to the total received signal power 
from source nodes. Note, however, that the present problem is not equivalent to the conventional extended 
network framework [10] due to quite different channel characteristics, and the main result is shown here 
in a somewhat different way by providing a simpler derivation than that of [10]. 

As illustrated in Fig. [T] let and Di denote the sets of sources and destinations, respectively, for 
the cut L in an extended network. We then take into account an approach based on the amount of power 
transferred across the network according to different operating regimes, i.e., path-loss attenuation regimes. 
As pointed out in [10], the information transfer from Sl to Dl is highly power- limited since all the nodes 
in the set Dl are ill-connected to the left-half network in terms of power. This implies that the information 
transfer is bounded by the total received power transfer, rather than the cardinality of the set D^. For 
the cut L, the total throughput T(n) for sources on the left is bounded by the (ergodic) capacity of the 
MIMO channel between Sl and Dl under time-varying channel assumption, and thus is given by 



where is the matrix with entries [11^]^^ = hki for i E SL,k e Dl, and Ql G C®^")''®^") is the 
positive semi-definite input signal covariance matrix whose fc-th diagonal element satisfies [QLjfcfe < P 
for k G Sl- 

The relationship (|9]) will be further specified in Theorem [T] Before that, we first apply the techniques 
of [26], [34] to obtain the total power transfer of the set Dl- These techniques involve the design of the 

^In [32], instead of simply taking the far-field assumption, the physical limit of wireless radio networks has been studied under certain 
conditions on scattering elements. Further investigation is also required to see whether this assumption is valid for underwater networks of 
unit node density in the limit of large number n of nodes. 



III. Cut-set Upper Bound 




(9) 
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optimal input signal covariance matrix Ql in terms of maximizing the upper bound (|9]) on the capacity. 
If the matrix has independent entries, each hki of which is a proper complex random variable [35], 
and has the same distribution as —hki for i £ SL,k G Di, then the optimal is diagonal, i.e., the 
maximum in ^ is attained with [QLjfcfc = P for k E Sl, where Ql is the diagonal matrix. We start from 
the following lemma. 

Lemma 1: Each element h^i of the channel matrix is a proper complex random variable, where 

ieSL,ke Dl. 

The proof of this lemma is presented in Appendix |A] It is readily proved that h^i has the same 
distribution as —h^i for all i and k since the random phases 9ki are uniformly distributed over [0,27r). 
Thus, using the result of Lemma [U we obtain the following result. 

Lemma 2: The optimal input signal covariance matrix that maximizes the upper bound (|9]) is unique 
and is given by the diagonal with entries [QiJ/cfc = P for k E Sl. 

We refer to Section III of [34] for the detailed proof. From Lemma |2l the expression ^ is then rewritten 

as 



Tin) < E 



< E 





P 



N if Air,. J) 
P 



Air,,J)Nif) 



(10) 



where the second inequality is obtained by applying generalized Hadamard's inequality [36] as in [10], 
[26]. The last two steps come from ([T]) and the fact that log(l + x) < x for any x, which is only tight as 
X is small. Note that the right-hand side (RHS) of (flOl) represents the total amount of received signal-to- 
noise ratio (SNR) from the set Sl of sources to the set Dl of destinations for a given cut L. To further 
compute ([TOl ). we define the following parameter 

P 



P. 



ki 



(11) 



for some constant Cq > independent of n, which corresponds to the total power received from the signal 
sent by all the sources i E Sl node k on the right (see ([T) and ([8])). For convenience, we now index 
the node positions such that the source and destination nodes under the cut L are located at positions 
i—ix + 1, iv) ^iid {^x^ ky), respectively, for i^, k^ = 1, ■ ■ ■ , y/n/2 and iy, ky = 1, ■ ■ ■ , y/n. The scaling 



resuh of ^ defined in (El) can then be derived as follows. 



Lemma 3: In an extended network, the term P^"* 



in (fTTI) is given by 

Pf = O [kl-aif)-'^) , 



(12) 



v^/2. 



where k^ represents the horizontal coordinate of node k E Dl for /c^. = 1, ■ 

The proof of this lemma is presented in Appendix |B] We are now ready to show the cut-set upper 
bound in extended networks. 
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Theorem 1: For an underwater regular network of unit node density, the total throughput T[n) is 
upper-bounded by 



where C2 > is some constant independent of n. 

Proof: From ([T]) and (fT0l)-(fT2]). we obtain the following upper bound on the total throughput T{n): 

^ V V p^*^^ 

N{f) ^ ^ ^ 

C3P^/n 1 



< 



< 



< 



< 



< 



iV(/) a(/) - 1 



where C3 and C4 are some positive constants independent of n, which is equal to (fTSl) . This completes the 
proof of the theorem. ■ 
We remark that this upper bound is expressed as a function of the absorption a{f) and the noise PSD 
N{f), whereas an upper bound for wireless radio networks depends only on the constant value a [10]. 
In addition, using ([3]), ([5]), and © in (fT3l) results in two other expressions on the total throughput 



r(„) . o I -^"-y;/')"^^ I (14, 



and 

for some positive constants ci and as shown in ([3]) and dH), respectively. Hence, from (fT4l) . it is seen 
that the upper bound is inversely proportional to the attenuation parameter a{f) and decays fast with 
increasing a(/), thereby leading to a highly power-limited network irrespective of the parameter a(/). 



B. Dense Networks 

In a dense network, it is necessary to narrow down the class of S-D pairs according to their Euclidean 
distance to obtain a tight upper bound. In this subsection, we derive a new upper bound based on hybrid 
approaches that consider either the sum of the capacities of the multiple-input single-output (MISO) 
channel between transmitters and each receiver or the amount of power transferred across the network 
according to operating regimes, similarly as in [27]. 

For the cut L, the total throughput T(n) for sources on the left half is bounded by the capacity of the 
MIMO channel between Sl and Dl, corresponding to the sets of sources and destinations, respectively, 
and thus is given by In the extended network framework, upper bounding the capacity by the total 
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received SNR yields a tight bound due to poor power connections for all the operating regimes. In a 
dense network, however, we may have arbitrarily high received SNR for nodes in the set that are 
located close to the cut, or even for all the nodes, depending on the path-loss attenuation regimes, and thus 
need the separation between destination nodes that are well- and ill-connected to the left-half network in 
terms of power. More precisely, the set Dl of destinations is partitioned into two groups Dl i and Dl^2 
according to their location, as illustrated in Fig. |2] Then, since Lemmas \T\ and [2] also hold for the dense 
network, by applying generalized Hadamard's inequality [36], we have 



T(n) < maxE 

Ql>0 



logdet ( In/2 + ^^HiQiHf 



NifY 



< E 

< E 
+E 



log det I In/2 + 



p 



Nif) 

logdet (l\n,,,\ + ^H,.,iHf,i^ 



logdet Iin, ,1 + 



\Dl,2\ ^y^Hi,2Hf 2 



(15) 



where H^, ^ is the matrix with entries [H^ = h^i for i E Sl, k E Dl^u and / = 1, 2. Note that the first 
and second terms in the RHS of (fT5l) represent the information transfer from Sl to D^ i and from 5*^ 
to Dl 2, respectively. Here, Dl i denotes the set of destinations located on the rectangular slab of width 
xlI \fn immediately to the right of the centerline (cut), where G {0, 1, ■ ■ ■ , -\/n/2}. The set Dl/i is 
given by Dl \ Dl^i- It then follows that |-Dl,iI = ^L\fn and \Dl/i\ = {y/n/2 — XL)\fri. 

Let T/(n) denote the /-th term in the RHS of ([HI) for / G {1, 2}. It is then reasonable to bound Txin) 
by the cardinality of the set Dl.\ rather than the total received SNR. In contrast, using the power transfer 
argument for the term T2(n), as in the extended network case, will lead to a tight upper bound. It is thus 
important to set the parameter xl according to the attenuation parameter a(/), based on the selection rule 
for Xl [27], so that only Dl,\ contains the destination nodes with high received SNR. To be specific, we 
need to decide whether the SNR received by a destination k E Dl from the set Sl of sources, denoted 
by P^^^/N{f), is larger than one, because degrees-of-freedom (DoFs) (also known as capacity pre-log 
factor) of the MISO channel are limited to one. If destination node k has the total received SNR greater 
than one, i.e., P^^'^ = u{N{f)), then it belongs to Dl,i. Otherwise, it follows that k E Dl,2- 



For analytical tractability, suppose that 



a(/) = e((l + eo)"') for/3G [0,cx)) 



(16) 



where eo > is an arbitrarily small constant, which represents all the operating regimes with varying 
(3. As before, let us index the node positions such that the source and destination nodes are located at 



iM+i and 



positions 

then obtain the following scaling results for P}"^ as shown below. 



respectively, for i^, k^ 

(k) 



, \/n/2 and iy, ky 



n. We 



■ L 

Lemma 4: In a dense network, the term P)^ ' in (fTTI) is given by 



P, 



(k) 




if 1 < a < 2 and fc^ = o {n^^^~^+') 
if a = 2 and k^ = o {n^/^-'^+') 



(17) 



and 



(k) 




max < 1 



(l+eo) 



,1/2-/3 



(l+^o) 



if k^ 
if kr 



(18) 
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for arbitrarily small positive constants e and eo, where k^/ y/n is the horizontal coordinate of node k G D12. 

The proof of this lemma is presented in Appendix O Although the upper and lower bounds for P|'^^ are 
not identical to each other, showing these scaling results is sufficient to make a decision on 0;^ according 
to the operating regimes. It is seen from Lemma |4]that when k^ = o (77,^/^"'^+''), P^^ does not depend on 
the parameter /3 (or a(/)), while for k^ = (n^/^~''+''), node k E 2 gets ill-connected to the left half 
in terms of power since P^'^^ decreases exponentially with n. More specifically, when k^ = o 
it follows that ^ = u{n''^) from (dl, resulting in ^ = uj{N{f)) due to N{f) = 0(1). In contrast, 
under the condition k^ = (n^/^"'^^^), it is observed from (fTTT ) that P^^ is exponentially decaying as 

a function of n, thus leading to P|'^^ = o{N{f)). As a consequence, using the result of Lemma |4l three 
different regimes are identified and the selected xl is specified accordingly: 

( if/3 = 

XL = < nV2-/3+. if < /3 < 1/2 (19) 
[0 if /3 > 1/2 

for an arbitrarily small e > 0. It is now possible to show the proposed cut-set upper bound in dense 
networks. 

Theorem 2: Consider an underwater regular network of unit area. Then, the upper bound on the total 
throughput T(n) is given by 

{0{nlogn) if (3 = 

O (n'-f^+nogn) ifO</3<l/2 (20) 

where e and eo are arbitrarily small positive constants, and 05 and /S are defined in dH) and (fT6l ), respectively. 

Proof: We first compute the first term Ti(n) in (fTSi) . focusing on the case for < /3 < 1/2 since 
otherwise T'i(n) = 0. Since the nodes in the set D^ i have good power connections to the left-half 
network and the information transfer to Di i is limited in bandwidth (but not power), the term Ti{n) 
is upper-bounded by the sum of the capacities of the MISO channels. More specifically, by generalized 
Hadamard's inequality [36], Ti{n) can be easily bounded by 

-"^^^°H'+a(/)Vy.iV(/)j 
<a;.v/^log^l + ^^J 

< c^XL^/nlogn (21) 

for some constant C5 > independent of n, where the last two steps are obtained from the fact that 
< a(/) < 1 and A^(/) tends to decrease polynomially with n from the relation in (|6]). The upper bound 
for the second term T2{n) in (fTSi) is now derived under the condition (3 E (0, 00). Similarly as in the 
steps of (flOl) . we have 



s EE 



Mni, f) 

P 



J2 (22) 
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which corresponds to the sum of the total received SNR from the left-half network to the destination set 
Dl 2- Hence, combining the two bounds (|2T]) and (l22l) along with the choices for specified in (fT9l) . 
we obtain the following upper bound on the total throughput T{n): 



Tin) < < 



if (3 = 
E ifO</3<l/2 

if /3 > 1/2 



< 



c^n log n 
1 V p(*^) 

N(f) ^ 
k&DL 

c^n log n 

c^n iog a -f- j^^j,^ l^k^=XL l^ky=i 

Ar(/) l^kx=l 2^ky=l 

c^n log n 
CqU^'^^" log n 

„(l+")/2 



1/2+1/2-/3 



1/2 



(l+^o) 



(l+£0)'= 



if /3 = 

if < /3 < 1/2 

if /3 > 1/2 



iV(/) 



Y^V"/2 1 



if /3 = 
if < /3 < 1/2 

if /3 > 1/2 



(23) 



for an arbitrarily small e > and some constant Cg > independent of n, where the equality comes from 



(fTTI) . The second inequality holds due to the fact that the term 



1/2 



fc^n/3-1/2 



tends to decay exponentially 



with n under the conditions < /3 < 1/2 and xl < < ^Jnl2, and thus the total is dominated by 
the information transfer T\{n) in (|2T)) . Now let us focus on the last line of (|23l) . which corresponds to 
the total amount of SNR received by all nodes for the condition /3 > 1/2. For /3 > 1/2, using the two 
relationships I© and (fT6] ) follows that 



y 



< 



^(l+«)/2 



^(l+")/2 

,(l+a+/3a5)/2 



< 



(3-1/2 



„/3-l/2 



u;(l) 



for some constant Cy > independent of n, where the second inequality holds due to (l + eo) 
under the condition. This coincides with the result shown in (|20l) , which completes the proof. ■ 
Note that the upper bound [10] for wireless radio networks of unit area is given by 0(nlogr2), which 
is the same as the case with /3 = (or equivalently a(/) = 0(1)) in the dense underwater network. Now 
let us discuss the fundamental limits of the network according to three different operating regimes shown 
in ([201). 

Remark 1: The upper bound on the total capacity scaling is illustrated in Fig. [3] versus the parameter 
/3 (logarithmic terms are omitted for convenience). We first address the regime /3 = (i.e., low path- 
loss attenuation regime), in which the upper bound on T{n) is active with x^ = \/n/2, or equivalently 
Dl i = Dl, while T2{n) = 0. In this case, the total throughput of the network is limited by the DoFs 
of the Q{n) x Q(n) MIMO channel between Sl and D^, and is roughly linear in the bandwidth, thus 
resulting in a bandwidth-limited network. In particular, our interest is in the operating regimes for which 
the network becomes power-limited as /3 > 0. In the second regime < /3 < 1/2 (i.e., medium path-loss 
attenuation regime), as pointed out in the proof of Theorem |2l the upper bound on T{n) is dominated by 
the information transfer from Sl to -Dl,i, that is, the term Ti(n) in (|2T]) contributes more than T2{n) in 
((22)) . The total throughput is thus limited by the DoFs of the MEMO channel between Sl and Dl,i, since 
more available bandwidth leads to an increment in Ti{n). As a consequence, in this regime, the network is 
both bandwidth- and power-limited. In the third regime (3 > 1/2 (i.e., high path-loss attenuation regime), 
the upper bound (|22l) is active with xl = 0, or equivalently Dl^2 = -Dl, while Ti{n) = 0. The information 
transfer to Dl is thus totally limited by the sum of the total received SNR from the left-half network to 
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the destination nodes in D^. In other words, in the third regime, the network is limited in power, but not 
in bandwidth. 

Note that the upper bound on T{n) decays polynomially with increasing (3 in the regime < /3 < 1/2, 
while it drops off exponentially when (3 > 1/2. In addition, two other expressions on the total throughput 
T{n) are summarized as follows. 

Remark 2: From (|6]) and (fT6l ). the upper bound and the corresponding operating regimes can also be 
presented below in terms of the attenuation parameter a(/): 

{0(n\ogn) ifa(/) = e(l) 

O (iSf') if «(/) = c^(l) and a(/) = O ((1 + e^)^) 

ifa(/)=c.((l + eo)v^). 

Note that as a(/) = w ((1 + cq)^), we also obtain 

^(l+a)/2 \ 



T(n) = O 



which is expressed as a function of the spreading factor a as well as the absorption a(/) and the noise 
PSD N(f). Using (|3]) and (|5]) further yields the following expression 

O(nlogn) if/ = 0(1) 
T(n) = l 0(^^^^^) if/ = c^(l)and/ = 0(ni/4) 



which represents the upper bound for the operating regimes identified by frequency scaling. 

IV. ACHIEVABILITY RESULT 

In this section, to show the order optimality in underwater networks, we analyze the achievable 
throughput scaling for both extended and dense networks with the existing transmission scheme, commonly 
used in wireless radio networks. Under an extended regular network, the conventional MH transmission is 
introduced and its optimal achievability result is shown. Under a dense regular network, we examine the 
operating regimes for which the achievable throughput matches the upper bound shown in Section IIII-BI 



A. Extended Networks 

The nearest-neighbor MH routing protocol [1] will be briefly described to show the order optimality. 
The basic procedure of the MH protocol under our extended regular network is as follows: 

• Divide the network into square routing cells, each of which has unit area. 

• Draw an line connecting a S-D pair. A source transmits a packet to its destination using the nodes 
in the adjacent cells passing through the line. 

• Use the full transmit power at each node, i.e., the transmit power P. 

The achievable rate of MH is now shown by quantifying the amount of interference. 

Lemma 5: Consider an extended regular network that uses the nearest- neighbor MH protocol. Then, 
the total interference power P/ from other simultaneously transmitting nodes, corresponding to the set 
/ C {1, ■ ■ ■ ,n}, is upper-bounded by 6(l/a(/)), where a{f) denotes the absorption coefficient greater 
than 1. 
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Proof: There are 8A; interfering routing cells, each of which includes one node, in the fc-th layer 1^ 
of the network as illustrated in Fig. |4] Then from ([1), (|7]), and ([8]), the total interference power P/ at each 
node from simultaneously transmitting nodes is upper-bounded by 



8P ^ 



< 



< 



y ' 

8P^ 1 
CO ^ a(/)^ 

cs 



«(/)' 

where Cq and Cg are some positive constants independent of n, which completes the proof. ■ 
Note that the received signal power no longer decays polynomially but rather exponentially with 
propagation distance in our network. This implies that the absorption term a(/) in ([T]) will play an 
important role in determining the performance. It is also seen that the upper bound on Pj does not depend 
on the spreading factor a. Using Lemma [51 it is now possible to simply obtain a lower bound on the 
capacity scaling in the network, and hence the following result presents the achievable rates under the 
MH protocol. 

Theorem 3: In an underwater regular network of unit node density. 



is achievable. 

Proof: Suppose that the nearest-neighbor MH protocol is used. To get a lower bound on the capacity 
scaling, the signal-to-interference-and-noise ratio (SINR) seen by receiver i E {1, ■ ■ ■ ,n} is computed as 
a function of the absorption a(/) and the PSD N(f) of noise rij. Since the Gaussian is the worst additive 
noise [37], [38], assuming it lower-bounds the throughput. Hence, by assuming full CSI at the receiver, 
from (HI), ([7]), and ([8]), the achievable throughput per S-D pair is lower-bounded by 

log(l + SINR) 

P/Mf)) 



> log 1 + 



> log 1 + 



N{f) + cs/a{f) 
C9P \ 

for some positive constants co, cg, and cg independent of n, where the second inequality is obtained 
from the relationship (|6]) between a(/) and N(f), resulting in N(f) = n (l/a(/)). Due to the fact that 
log(l + x) = (loge)2; + 0{x'^) for small x > 0, the rate of 

Q 



is thus provided for each S-D pair. Since the number of hops per S-D pair is given by 0{y/n), there 
exist n(y/n) source nodes that can be active simultaneously, and therefore the total throughput is finally 
given by (|24)) . which completes the proof of the theorem. ■ 
Now it is examined how the upper bound shown in Section IIII-AI is close to the achievable throughput 
scaling. 
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Remark 3: Based on Theorems \T\ and |3l it is easy to see that the achievable rate and the upper bound 
are of exactly the same order. MH is therefore order-optimal in regular networks with unit node density 
for all the attenuation regimes. 

We also remark that applying the hierarchical cooperation strategy [10] may not be helpful to improve 
the achievable throughput due to long-range MIMO transmissions, which severely degrade performance 
in highly power-limited networks]^ To be specific, at the top level of the hierarchy, the transmissions 
between two clusters having distance 0{^yn) become a bottleneck, and thus cause a significant throughput 
degradation. It is further seen that even with the random phase model, which may enable us to obtain 
enough DoF gain, the benefit of randomness cannot be exploited because of the power limitation. 

B. Dense Networks 

From the converse result in Section ITII-Bi it is seen that in dense networks, there exists either a bandwidth 
or power limitation, or both, according to the path-loss attenuation regimes. Based on the earlier studies [1], 
[10], [27], [39] for wireless radio networks, it follows that using MH routing is preferred at power-limited 
regimes, while the HC strategy may have a better performance at bandwidth-limited regimes. Thus, existing 
schemes need to be used carefully, depending on operating regimes. 

In this subsection, the nearest-neighbor MH routing in [1] is described with a slight modification. The 
basic procedure of the MH protocol under our dense regular network is similar to the extended network 
case, and is briefly described as follows: 

• Divide the network into n square routing cells, each of which has the same area. 

• Draw a line connecting an S-D pair. 

• At each node, use the transmit power of 



the received SNR at each node is bounded by 1 (note that having a higher power is unnecessary in terms 
of throughput improvement). 

The amount of interference is now quantified to show the achievable throughput based on MH. 

Lemma 6: Consider a dense regular network that uses the nearest-neighbor MH protocol. Then, the 
total interference power P/ from other simultaneously transmitting nodes, corresponding to the set / C 
{1, ■ ■ ■ , n}, is bounded by 



for an arbitrarily small eo > 0, where and (3 are defined in dD and (fT6l ). respectively. 

The proof of this lemma is presented in Appendix |Dl From (|6]) and ([T6l) . we note that when /3 = 1/2, 
it follows that P/ = 0{N{f)), i.e., P/ is upper-bounded by the PSD A^(/) of noise. Using Lemma [6l a 
lower bound on the capacity scaling can be derived, and hence the following result shows the achievable 
rates under the MH protocol in a dense network. 

'in wireless radio networks of unit node density, tlie liierarchical cooperation provides a near-optimal throughput scaling for the operating 
regimes 2 < a < 3, where a denotes the path-loss exponent that is greater than 2 [10]. Note that the analysis in [10] is valid under the 
assumption that a is kept at the same value on all levels of hierarchy. 




O 





(25) 
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Theorem 4: In an underwater regular network of unit area, 



Tin) 



^ ' |„(1/2-;3)(2-q)_ 
/ ^(l+a+;3a5)/2 ' 



if < /3 < 1/2 

if /3 = 1/2 
if /3 > 1/2 



(26) 



is achievable. 

Proof: Suppose that the nearest-neighbor MH protocol described above is used. Then, from ([T]), the 
received signal power from the desired transmitter is given by 



P min < 1 



a(/)V^jV(/) 1 
n"/2 / 



n 



a/2 



which can be rewritten as 



Co(l + Co) 

P 



n/3-1/2 



mm 



;i+eo) 



,/3-l/2 



' ^{Q+/3a5)/2 

if < /3 < 1/2 
if /3 > 1/2 



(27) 



with respect to the parameter (3 using ([61) and (fT6l ). A lower bound on the throughput is now obtained 
using (|25l) and (ITTI) . By assuming the worst case noise, which lower-bounds the transmission rate, and 
full CSI at the receiver, the achievable throughput per S-D pair is then lower-bounded by 

log(l + SINK) 



log 1 + 



Pr 



Nif) + Pi 



n log 1 + 



1 



1 + 



„(a + /3a5)/2 



/2-/3)(2-c«) 



,logn} 



if < /3 < 1/2 

if /3 = 1/2 
if /3 > 1/2 



1](1) 



{n(i/2-;3)(2-Q)^log„]. 



„(a+/3a5)/2 



(l+^o) 



„/3-l/2 



if < /3 < 1/2 

if /3 = 1/2 
if /3 > 1/2, 



where the second equality holds since A^(/) = 6(n"^"5/2) tj^yg ^ uj{N{f)) for < /3 < 1/2, 
P/ = e(P^) = e(iV(/)) for /3 = 1/2, and P/ = o(A^(/)) for /3 > 1/2. The last equality comes from the 
fact that log(l + x) = (loge)x + O(x^) for small x > 0. Since there are ri{y/n) S-D pairs that can be 
active simultaneously in the network, the total throughput is finally given by (l26l) . which completes the 
proof. ■ 
Note that the achievable throughput [1] for wireless radio networks of unit area using MH routing is 
given by fl{^/n), which is the same as the case for which /3 = 1/2 (or equivalently a(/) = 6 ((1 + eo)^)) 
in a dense underwater network. The lower bound on the total throughput T(n) is also shown in Fig. |3] 
according to the parameter (3. From this result, an interesting observation follows. To be specific, in the 
regime < /3 < 1/2, the lower bound on T(n) grows linearly with increasing (3, because the total 
interference power P/ in (|25l) tends to decrease as /3 increases. In this regime, note that P/ = ^l(Pr). 
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Meanwhile, when /3 > 1/2, the lower bound reduces rapidly due to the exponential path-loss attenuation 
in terms of increasing (3. 

In addition, similarly as in Section IIII-Bl two other expressions on the achievability result are now 
summarized as in the following. 

Remark 4: From ^ and (fT6] ). the lower bound on the throughput T{n) and the corresponding operating 
regimes can also be presented below in terms of the attenuation parameter a(/): 



T{n) 



^ ( f , ,,y^\; r- r I if = ^(1) and a(/) = o ((1 + eo)^) 

a(^i) if a(/) = e ((1 + e„)^) 



Furthermore, using (|3]) and (|5} follows that 



Tin) 



^ ( r r2^2/^ if / = ^(1) and / = o (n^/^) 

n (V^) if / = e (n^/^) 



which represents the lower bound for the operating regimes obtained from the relationship between the 
frequency / and the number n of nodes. 

Now let us turn to examining how the upper bound shown in Section IIII-BI is close to the achievable 
throughput scaling. 

Remark 5: Based on Theorems |2]and|4l it is seen that if /3 > 1/2, then the achievable rate of the MH 
protocol is close to the upper bound up to rf for an arbitrarily small e > (note that the two bounds are of 
exactly the same order especially for (3 > 1/2). The condition /3 > 1/2 corresponds to the high path-loss 
attenuation regime, and is equivalent to a(/) = i7 ((1 + cq)^) or / = ^7 (n^/'^). Therefore, the MH is 
order-optimal in regular networks of unit area under the aforementioned operating regimes, whereas in 
extended networks, using MH routing results in the order optimality for all the operating regimes. 

Finally, we remark that applying the HC strategy [10] does not guarantee the order optimality in 
the regime < /3 < 1/2 (i.e., low and medium path-loss attenuation regimes). The primary reason is 
specified under each operating regime: for the condition /3 = 0, following the steps similar to those of 
Lemma |6l it follows that Pj = uj{Pr) at all levels of the hierarchy, thereby resulting in SINR = o(l) 
for each transmission (the details are not shown in this paper). It is thus not possible to achieve a 
linear throughput scaling. Now let us focus on the case where < /3 < 1/2. At the top level of the 
hierarchy, the transmissions between two clusters having distance 0(1) becomes a bottleneck because of 
the exponential path-loss attenuation with propagation distance. Hence, the achievable throughput of the 
HC decays exponentially with respect to n, which is significantly lower than that in (|26] ). 



V. Extension to Random Networks 

In this section, we would like to mention a random network configuration, where n S-D pairs are 
uniformly and independently distributed on a square of unit node density (i.e., an extended random 
network). 

We first discuss an upper bound for extended random networks. A precise upper bound can be obtained 
using the binning argument of [10] (refer to Appendix V in [10] for the details). Consider the same cut L, 
which divides the network area into two halves, as in the regular network case. For analytical convenience, 
we can artificially assume the empty zone El, in which there are no nodes in the network, consisting of 
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a rectangular slab of width < c < -^jf^, independent of n, immediately to the right of the centerline 
(cut), as done in [27] (see Fig. [5])0 Let us state the following lemma. 

Lemma 7: Assume a two dimensional extended network where n nodes are uniformly distributed. When 
the network area is divided into n squares of unit area, there are fewer than logn nodes in each square 
with high probability. 

Since the result in Lemma |7] depends on the node distribution but not the channel characteristics, the 
proof essentially follows that presented in [4]. By Lemma |7l we now take into account the network 
transformation resulting in a regular network with at most log n and 2 log n nodes, on the left and right, 
respectively, at each square vertex except for the empty zone (see Fig. |5]). Then, the nodes in each square 
are moved together onto one vertex of the corresponding square. More specisely, under the cut L, the node 
displacement is performed in the sense of decreasing the Euclidean distance between source node i £ Sl 
and the corresponding destination k G D^, as shown in Fig. |5J which will provide an upper bound on P^^ 
in (fTTI) . It is obviously seen that the amount of power transfer under the transformed regular network is 
greater than that under another regular network with at most log n nodes at each vertex, located at integer 
lattice positions in a square region of area n. Hence, the upper bound for random networks is boosted by 
at least a logarithmic factor of n compared to that of regular networks discussed in Section lUl 

Now we turn our attention to showing an achievable throughput for extended random networks. In 
this case, the nearest- neighbor MH protocol [1] can also be utilized since our network is highly power- 
limited. Then, the area of each routing cell needs to scale with 2 log n to guarantee at least one node 
in a cell [1], [6]0 Each routing cell operates based on 9-time division multiple access to avoid causing 
large interference to its neighboring cells [1], [6]. For the MH routing, since per-hop distance is given by 
( A/log n), the received signal power from the intended transmitter and the SINR seen by any receiver 
are expressed as 

CiqP 

(logn)°/2a(/)'^v'i°g" 

and 



(logn)"/2a(/)'5v^A^(/) 



respectively, for some co nstants ci p > and 5 > \ f2 indep endent of n. Since the number of hops per 
S-D pair is given by 0{^Jn/ \ogn), there exist U{^JnJ\ogn) simultaneously active sources, and thus the 
total achievable throughput T{n) is finally given by 



Tin) = Q 



n 



1/2 



(log n)("+i)/2a(/)'5vT^A^(/) 

for some constant 5 > \/2 independent of n (note that this relies on the fact that log(l + x) can be 
approximated by x for small x > 0). Hence, using the MH protocol results in at least a polynomial 
decrease in the throughput compared to the regular network case shown in Section |IV0 This comes from 
the fact that the received signal power tends to be mainly limited due to exponential attenuation with 
transmission distance ©(Vlogn). Note that in underwater networks, randomness on the node distribution 
causes a huge performance degradation on the throughput scaling. Therefore, we may conclude that the 

'"Although this assumption does not hold in our random configuration, it is shown in [27] that there exists a vertical cut such that there 
are no nodes located closer than < c < ^^1/4 on both sides of this cut when we allow a cut that is not necessarily linear. Such an 
existence is proved by using percolation theory [4], [40]. This result can be directly applied to our network model since it only relies on 
the node distribution but not the channel characteristics. Hence, removing the assumption does not cause any change in performance. 

"when methods from percolation theory are applied to our random network [4], [40], the routing area constructed during the highway 
phase is a certain positive constant that is less than 1 and independent of n. The distance in the draining and delivery phases, corresponding to 
the first and last hops of a packet transmission, respectively, is nevertheless given by some constant times log n, thereby limiting performance, 
especially for the condition a{f) = a;(l). Hence, using the routing protocol in [4] indeed does not perform better than the conventional MH 
case [1] in random networks. 

'"In terrestrial radio channels, there is a logarithmic gap in the achievable scaling laws between regular and random networks [1], [25]. 
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existing MH scheme does not satisfy the order optimality under extended random networks regardless of 
the attenuation parameter a{f). 



VI. Conclusion 

The attenuation parameter and the capacity scaling laws have been characterized in a narrow-band 
channel of underwater acoustic networks. Provided that the frequency / scales relative to the number n of 
nodes, the information-theoretic upper bounds and the achievable throughputs were obtained as functions 
of the attenuation parameter a(/) in regular networks. In extended networks, based on the power transfer 
argument, the upper bound was shown to decrease in inverse proportion to a{f). In dense networks, the 
upper bound was derived characterizing three different operating regimes, in which there exists either a 
bandwidth or power limitation, or both. In addition, to show the achievability result, the nearest-neighbor 
MH protocol was introduced with a simple modification, and its throughput scaling was analyzed. We 
proved that the MH protocol is order-optimal in all operating regimes of extended networks and in power- 
limited regimes (i.e., the case where the frequency / scales faster than or as n^/'') of dense networks. 
Therefore, it turned out that there exists a right frequency scaling that makes our scaling results for 
underwater acoustic networks to break free from scaling limitations related to the channel characteristics 
that were described in [23]. Our scaling results were also extended to the random network scenario, where 
it was shown that the conventional MH scheme does not satisfy the order optimality for all the operating 
regimes. 



Appendix 

A. Proof of Lemma [7] 

The following definition is used to simply provide the proof. 

Definition 1 [35]: A complex random variable Y is said to be proper if Sy = 0, where Sy, called the 
pseudo-covariance, is given by E[(Y — E[Y]y]. 

Since the {k,i)-th element of the channel matrix Hj^ is given by ([8]), it follows that 

1 



E[{hki-E[hk^]r] 



-E 



From the fact that 

E [e^'"^] = E [cos(0,,) + J sin(4,)] = 
due to uniformly distributed O^i over [0, 2n], we thus have 

.2. 1 



E[{hk^-E[hk.]y] 



1 

0, 



E[cos{2eki)+jsm{2eki)] 



(28) 



which complete the proof, because (|28l) holds for all i e Sl and k E Dl- 



B. Proof of Lemma \3\ 

An upper bound on P^^ can be found by using the node-indexing and layering techniques similar to 
those shown in Section VI of [39]. As illustrated in Fig. |6l layers are introduced, where the z-th layer l[ of 
the network represents the ring with width 1 drawn based on a destination node k E Dl, whose coordinate 
is given by {k^, ky), where i G {1, ■ ■ ■ , y/n). More precisely, the ring is enclosed by the circumferences 
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of two circles, each of which has radius k^ + i and k^ + i — 1, respectively, at its same center (see Fig. |6]). 
Then from (fTTT ). the term P^'^ is given by 



.w t^i ((t + 1 - 1)= + (<y - k,r-r" a(/)V('-+'--i)"+<'.-'.)' 

It is further assumed that all the nodes in each layer are moved onto the innermost boundary of the 
corresponding ring, which provides an upper bound for P^\ From the fact that there exist &{kj. + i) 
nodes in the layer /• since the area of /• is given by TT{2kx + 2i — 1), p|^^ is then upper-bounded by 



p{k) < :^ cii(z^ + l) 
2cuP ^ 1 

2^ n(fV 



^ 2ci,p ( 1 r 1 

^ C12P 



for some positive constants cq, ch, and C12 independent of n, where the fourth inequality holds since 
a{f) > 1, which finally yields (fT2l) . This completes the proof. 

C. Proof of Lemma |?] 

Upper and lower bounds on P^''^ in a dense network are derived by basically following the same node 
indexing and layering techniques as those in Appendix |B] We refer to Fig. |6] for the detailed description 

(note that the destination nodes are, however, located at positions ^ j in dense networks). Similarly 
to the extended network case, from (fTTT ). the term P^^^ is then given by 



Co ,fr; ((i, + A. - 1)2 + (i, - 

First, focus on how to obtain an upper bound for P^''^ . Assuming that all the nodes in each layer are 
moved onto the innermost boundary of the corresponding ring, we then have 

p(fc) ^ P^ ^ + IJ 

^ - Co ^'°«(/)"^ 



y I 



= Ci2Pn"/2 \^ ^ , , (29) 

for some positive constants cq, cn, and C12 independent of n and an arbitrarily small eo > 0, where 
the equality comes from the relationship (fT6] ) between a{f) and (3. We first consider the case where 
kx = o (n^/^^^^'') for an arbitrarily small e > 0. Under this condition, from the fact that the term 
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in the RHS of (|29l) is dominant in terms of upper-bounding P^'^'' for i' = k^, ■ ■ ■ , y/n, (|29l ) is further 
bounded by 



1 rV^ 1 

< Ci2Pn°/2 + / dx 

— ' Ua—l I rra — 1 



which yields P^*^^ = O (ri"/2(y^)2-") = for 1 < a < 2 and Pf ^ = O (rilogn) for a = 2. When 



/cj; = f2(n^/^ ^'*"^), the upper bound (1291 ) for P^^^ is dominated by the term (1 + eo)*"^ ^ , and thus is 
given by 



Pf) <Ci2Pn"/2y i ^ 



1 /"^ 1 

^ ci2Pn"/2 ^ :7-^ + / —^dx 



.1 + Co) (1 + eo) 

1 rVn/k^ ^ 



J- + CoJ U + eoj 

<-^^^^^!^--{l>^^/'-'} (30) 
for some constant cis > independent of n, which is the last result in ([TT] ). 

(k) 

Next, let us turn to deriving a lower bound for P^ . Since each layer has at least one node that is onto 
the innermost boundary of the corresponding ring, the lower bound similarly follows 

p„a/2 fe-+v^/2-l . 

pik) ^ Pn^ ^ 1 



ci2Pn"/2 ^-^JZTH- (31) 

i'=fc, ^''^ (1 + eo) 



For the condition = o^n^^"^ ^~^''), (l3T1) is represented as 



Pf ) > Ci2Pn-/2 ^ 



A:;c+V"/2-l ^ 



i =kx 



(1 + eo) 



2fc^-l 

> Ci2Pr2°/2 ^ 



^ Ci2Pn°/2 



(1 + eo) 



y - 



i =hx 



~ n^'k^-'^ 
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for an arbitrarily small e' > and some constant C14 > independent of n, where the third inequality 



holds due to (1 + eo 



0{rf'). On the other hand, similarly as in the steps of (|30] ). the condition 



kx = '^^^) yields the following lower bound for P^^^: 



> 



>Ci2P 
>Ci5P 
> - 



i'=k^ (1 + Cq 

1 



+ 



:i + eo) " Jk,+i 



1 



;i + eo 



>xn/3-i/2 



dx 



N h n/3-1/2 



+ 



1/2-/3 



:i+eo 



.(fc^+lK-i/a 



max < 1, 



n 



1/2-/3 



1/3-1/2 



+ I (1 + eo) 

some constant C15 > independent of n, which finally complete the proof of the lemma. 

D. Proof of Lemma |^ 

The layering technique illustrated in Fig. |4]is applied as in the extended network case. From ([T), the 
total interference power P/ at each node from simultaneously transmitting nodes is then upper-bounded 
by 



k=l 



8PW2min{l,^(^l^^} 



Co 



y I . 



(32) 



Using ^ and (fT6l) . the upper bound (|32] ) on P/ can be expressed as 



Pi < CiePn^/^min 



eoj 



,,,8-1/2 



^(a+/3a5)/2 



E 



< 



C16-P y^V" 1 

?1^^572 Z^fc = l jr.a-l(]^_|_^|j)fen''-l 

cief sr^Vn 1 

„/3a5/2 Z^fc = l (l+eo)': 



72" 



< 



CiePn"/^ 1 

„/3a5/2 Z^fc=l 
n'3°5/2 (l+eo)-l 

CiePn"/^ 1 



1 



fcn/3-1/2 



(l+^o) 

1 

/c<»-l(l+eo)*"''"^''^ 



„/3-l/2_ 



< 



„/3a5/2 Z^fc = l 



fc°=-l(l+eo)'="''^^''^ 



l/3<'5/2 



x/3-1/2 



.=1 fc°-Hi + eo 

if < /3 < 1/2 
if /3 = 1/2 
if /3 > 1/2 

if < /3 < 1/2 
if /3 = 1/2 
if /3 > 1/2 

if < /3 < 1/2 
if /3 = 1/2 
if /3 > 1/2 



ifcn/3-i/2 
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for some positive constants Cig and Cn independent of n. Based on the argument in Appendix O when 
< /3 < 1/2, it follows that 



>^ 1 1 ^ 1 



fc=l fe=ni/2-/3 



ni/2-,9 



< I 1 + / -^dx 
n(V2-«(-i) [yTTo + y„v2-. (1 + 60)-"-^^^'^"' 

r.ni/2-/3 ^ 2 ^1/2-/3 



< 



4^(l/2-/3)(2-a) ifl<c,<2 

log n if a = 2, 



„.ax|n(i/2-'3)(2~"),logn| ~ 

which results in Pj = O [ — — ^ ^^^^72 — ] ■ This completes the proof of the lemma. 
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Fig. 1. The cut L in a two-dimensional extended regular network. Sl and Dl represent the sets of source and destination nodes, respectively. 
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Fig. 2. The cut L in a two-dimensional dense regular network. Sl and Dl represent the sets of source and destination nodes, respectively, 
where Dl is partitioned into two groups Dl.i and Dl,2- 
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Fig. 3. Upper (solid) and lower (dashed) bounds on the capacity scaling T{n). 
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Fig. 4. Grouping of interference routing cells in extended networks. The first layer l\ represents the outer 8 shaded cells. 
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Fig. 5. The node displacement to square vertices, indicated by arrows. The empty zone with width constant c is assumed for simphcity. 
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Fig. 6. Grouping of source nodes in extended networks. There exist 0(fca;) nodes in the first layer V^. This figure indicates the case where 
one destination is located at the position (kx, ky). The source nodes are regularly placed at spacing 1 on the left half of the cut L. 



